function [phi, C]=reg(x,p)

X=[];

for j=1:p

    X=[X x(:,p-j+1:end-j)'];

end

y=x(:,p+1:end)';

phi=(X'*X)\X'*y;

res=y-X*phi;
C=cov(res);
